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ABSTRACT 

This paper presents the theory and simulation of the derivative line shape function of the EPR orthorhombic 
Jahn-Teller spectra of Cu 2+ in Cd 2 (NH 4 )2(S04)3 single crystals. Computer simulation of the Cu 2+ EPR spectra 
using the spectrum fitting technique has been carried out in the temperature range of 15-180 K. The results 
show excellent fit between the simulated and observed spectra. The insufficient resolution of the spectra at 15 K 
is attributed to the coupling of copper isotopes and its orthorhombicity to phonon-induced reorientation rate or 
quantum tunnelling of Cu 2+ ions among the three JT potential wells. The splitting processes in the 3d 9 orbital 
states ofCu 2+ has been explained and in conjunction with the demands of the criterion of pure static JT systems, 
the random strain splitting is found to be 0.0176 cm' 1 and 0.0342 cm' 1 at 15 K and 180 K respectively. 
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I. INTRODUCTION 

Jona and Pepinsky (1956) discovered the Langbeinite family of compounds with general chemical 
formula (X + ) 2 (Y 2+ ) 2 (S0 4 )3 where X + is ammonia or monovalent metal and Y is a divalent metal. According to 
Babu et al (1984), many authors have carried out a number of investigations on the microscopic properties of the 
langbeinite family of crystals. The compound cadmium ammonium sulphate, Cd 2 (NH4)2(S0 4 )3 abbreviated as 
CAS is isomorphous to potassium magnesium sulphate, K 2 Mg2(S0 4 )3 abbreviated as PMS and both are 
members of the langbeinite family of crystals. The detailed x-ray data regarding the atomic positions in CAS 
does not seem to be available in the literature but that of PMS have been reported by Zemann and Zemann 
(Babu et al, 1984). Also they stated that some authors have shown that the langbeinite family of compounds in 
cubic space group P 2 /3 with a=10.35 A have four molecules per unit cell. The structure consists of a group of 
(S0 4 ) 2 " tetrahedral and Cd 2+ metal ions. There are two each crystallographically non -equivalent Cd 2+ and (NH4) + 
sites. Each Cd 2+ is surrounded by six oxygen atoms which form a slightly distorted octahedron. Yerima (2005, 
2007) stated that EPR studies of Misra and Korezak in 1986 using Mn 2+ as a probe revealed a phase transition at 
94.5 K whose mechanism was attributed to the freezing out of the rotation of the (S0 4 ) 2 " ion. Also he stated that 
EPR studies of Mouli and Sastry in 1962 using Cu 2+ probe at room temperature and 77 K yielded eight poorly 
resolved hyperfine lines in a general direction and a set of four unresolved hyperfine lines in any 
crystallographic plane. Their computed values of the g-factor showed that g//>g± and that the A max is falling 
along gmin. They concluded that Cu 2+ ions in this system may be in a compressed octahedral position or entered 
into the system inter stitially rather than substitutionally. The shortcomings of the study of Mouli and Sastry 
include its limitation to 77 K and its failure to focus on JT effect that could be associated with Cu 2+ in this 
crystal as earlier hinted by Babu et al in 1984. Raman spectroscopy performed by Rabkin et al (1981) provided 
substantial information on the structural phase transition in this crystal. They observed that CAS exhibits a 
structural phase transition on lowering the temperature from space group p 2 /3 to p 2 i at about 95 K, the lower the 
temperature phase being ferroelectric. Babu et al (1984) in their summary stated that Bhat et al (1973) EPR 
spectra of Mn 2+ in the two phase groups of CAS observed at room temperature and at liquid nitrogen were 
slightly different, particularly in the magnitudes of their zero-field splitting, as a result of small orthorhombic 
component in the low temperature phase. The differences in the spectra indicated a phase transition from high 
temperature phase p 2/3 to a low temperature phase p 2i . In another outlook, Yerima (2005) stated that Ng and 
Calvo in 1975 observed in their EPR studies using Mn 2+ in the temperature range 300-77 K observed definite 
change in the spectral pattern from that of room temperature to that of liquid nitrogen temperature due to a phase 
transition from p 2/3 to a space group of lower symmetry at low temperature. They did not analyze the spectra nor 
determine the phase transition temperature due to the complexity of the spectra at liquid nitrogen temperature. 
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The study of phase transition in CAS using VO 2 " as a probe in the temperature range of 573-77 K 
showed a complex spectrum, in which Babu et al (1984) were unable to identify the phase transition temperature 
and they suggested that the lower temperature phase stabilizes at temperatures higher than 95 K. This was due to 
a local stabilization of the low symmetry phase by vandyl ions. They also found that the data of their experiment 
contained Cu 2+ ion as contaminant and thus concluded that the complication in their results could have been due 
to a possible JT effect associated with Cu 2+ ion in this crystal. Unfortunately, they did not carry out further 
investigation to establish or rule out the possibility of JT effect due to Cu 2+ ion in this system. 

De (2010) said in order to ascertain JT effect of Cu 2+ in CAS Oguama carried out detailed EPR study 
of Cu 2+ doped in CAS at various temperatures in the range of 300-15 K with angular variations usually at steps 
of 5° and 0.5° near JT extrema points in three mutually perpendicular planes of CAS single crystals. The 
Cu 2+ :CAS spectra consists of finely resolved hyperfine lines in the low magnetic field region at all temperatures 
except 15 K while those in the high magnetic field region are poorly resolved at high temperatures. The complex 
or insufficient resolved splitting patterns of the spectra constitute a major problem in the extraction of the 
relevant spectroscopic parameters (resonance frequency, coupling constants, linewidths) from the spectra as 
starting points in the computer simulation of the spectra. Therefore, a straightforward analysis is not possible in 
general (neither for the human expert nor for a conceivable computer expert system (kirste, 1992)). It is against 
this background, several methods (spectrum contraction or elimination hyperfine splitting constants, extraction 
of hyperfine splitting constants, correlation method, maximum method, significance plots and "roll-up" 
transformation) have been employed in the analysis of high resolution EPR spectra (Kirste, 1992; Wu, 2006). 
These methods may be more or less useful in providing estimates for the hyperfine coupling constants. On the 
other hand, the set of relevant spectroscopic parameters extracted from the spectra can be verified by spectrum 
simulation and comparison with the experimental spectrum which is the method adapted in this paper. We have 
adapted the spectrum fitting because it is considered an indispensable component in any scheme automated 
spectrum analysis. 

Theory and Method of Computer Simulation of the Line Shape Function of Cu 2+ in CAS 

The detailed procedure of the experiment for determining Cu 2+ spectra in CAS is available in Oguama 
(1997) and De (2010). The Oguama (1997) spectra of Cu 2+ in this system at various temperatures were acquired 
for the computer simulation. The expression for the line shape function required for the computer simulation of 
the EPR spectra is derived from the spin Hamiltonian as follows: Thus, if the magnetic field is oriented along z- 
direction (100) of one of the orthorhombic JT distortions of the Cu 2+ in CAS, the field H z sees the g z as well as 
g x , g y tensor axes simultaneously unlike non-JT ions, and the spin Hamiltonian H s becomes for H along z, i.e. 
H x =H y =0. 

H s = ${H s S sSz + H s S xSx + H z S y3y ) + A x (S S I Z + S X I X +■ S y I y ) + A x (2S X I X -f V*> + Vy + 

1 

The expectation values of the Zeeman interaction for Cu 2+ can be deduced taking into consideration the tunnel 
splitting equations of d 9 ions given by (Boatner and Reynolds, 1977) 

ff 11 = r + p fl ,&s,+ji l s,4 

*22 = r + ^^ + ^v* 2 

H 22 = T + pg ± g z Ss + A ± S S I S 
where the terms with subscripts x or y in equation (1) vanishes i.e. H x =Hy=0, etc. The effective Hamiltonian for 
Cu 2+ is obtained by adding equations (2) bearing in mind the orthorhombic g and A tensors leading to 

H s =3T+ fi{HAfh + BA9x + HsSytfy) + ^ fr*4 + Vi + Vy) 

+ S S Q + A,{2S y I y + VJ 3 
where 3r is the tunnel splitting. The knowledge of the ratio of the random strain splitting, 5, to the tunnel 
splitting, 3r, tells us the type of JT effect present in the system in question. 

Considering the hyperfine interactions as perturbations up to second-order instead of a single transition between 
the non-JT Kramer's doublet, we have three transitions for magnetic field parallel z direction (in the case of a 
single Cu 2+ per unit cell in the static JT regime) for JT Kramer's doublet as follows (assuming that the principal 
axis of the g-tensor ellipsoid coincides, in orientation, with the principal axes of the hyperfine terms): 
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Jiv = faH„ + A Z M, + (JQ + 1] - Af/ ] 



(J(J + 1)-MP 



hv = pg y R zy +A y M } + V Q + l) - JM?) 



4 
5 
6 



Now considering only the first order terms in equations 4-6 the function describing line shapes in a non-JT 
complex for H along any g direction may be given by V(H} = H — + AM 3 . Thus, following the theory of the 

orthorhombic JT complex and considering transitions in equations 4-6, the resultant line shape function will be 
simultaneously governed by three functions given by 



P3y 



Thus the expression for the derivative EPR line shape can be written as 



J + 71 71 + 



+fl 3 } J 



10 



where fj. = g% (fflj - g^)cas z fp + g* + 

P* = Slisi ~ 8$™s*<p + g y + P 3 = glXgi - fl£ W* + + 
Jf j_ = ~ &)t ± ' = ~^t^-^2 = ~ ^t^G 2 = gxcas z <p + gy(l — x 2 '), x = cos9 = 1, = g x cos<p f 6 2 = 

with cot(=i,2,3) being the peak-to-peak linewidth (AH PP ). Finally integrating equation (10), the derivative line 
shape function of EPR spectra of Cu 2+ manifesting static JT effect with orthorhombic g and A tensors in CAS is 
given by 



{2a / : } 



11 



where f L = gl - ^ + ng$ t f 2 = g 2 (g* - g 2 ) ^ + nr^* f 3 = g 2 (g* - + ir fl *, a is the fractional 

relative abundance of copper isotopes a(Cu 63 ) = 0.69 and a(Cu 65 ) = 0.31 ; g x , g y , g z are g-factors along the x-, y- 
, and z-axes respectively; H is the applied magnetic field; A x , A y , A z are the hyperfine constants in the x-, y-, 
and z-axes respectively; Ri, R 2 , R3 are the linewidth parameters in the x-, y-, and z-axes respectively; Mi=-3/2, - 
1/2, 1/2, 3/2 and P is the Bohr magneton. The terms in the curly brackets represent the magnification factor of 
the line shape while the terms in square brackets determine the shape of the EPR line. The spectra of Cu 2+ in 
CAS at low temperatures consists of three sets of four resolved hyperfine lines in both the g// and g± directions. 
This implies that at the point of inflection the derivative line shape function is zero or cuts the H-axis four times 
in each of the three sets of four resolved lines. Thus the values of H corresponding to the points of inflection or 
zero F(H) can be determined. This means that the terms in square brackets must be zero i.e. H — — + AM = 0. 

Using two values of H(Hi and H 2 ) and M 2 (m 7i cmri Mj 2 ) corresponding to inflection points where F(H)=0, the 

corresponding expression for hyperfine constants is given by ^[[[=1^33 = ^* H± — ' 1 where n is the 



il-L 



number of complete resolved hyperfine lines and n-1 gives the number of hyperfine constants. For example, for 
Cu 63 M h = -3/2 and M h = 3/2 so that n-1 = M h - M h = 3 which is equivalent to n = 4. 
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On the other hand, generally the g values are obtained by noting the middle field point H^f of the 

hv 

hyperfine lines corresponding to M I± = —1/2 and M; 2 = 1/2 and using the relation &iti=x J y J s') = JZ — • The 

P**TTlf 

calculated A and g values and estimated value of peak-to-peak linewidths AH pp .(i=x,y,z) were used as starting 

points in the computer simulation process. These values of g, A and AH pp were substituted in the expression of 

F(H) and the resulting line shape was compared or matched with observed EPR line. The values of g, A and 
Atfpp were varied by small amounts in the ranges ±0.005-0.01, ±0.5-1 G, and ±5-10 G respectively around the 

initial values until the line that best fits the observed EPR line was obtained. The presence of M-shape at the 
first crest and W-shape at the last trough corresponding to g// lines are assumed to be due to coupling between 
Cu 63 and Cu 65 isotopes at low temperatures (Yerima, 2007). Now, if Cu 63 was the only one present in the sample 
of CAS we would expect only one set of four finely resolved hyperfine g /f lines corresponding to Mi=-3/2, -1/2, 
1/2, 3/2 as opposed to six lines for only Cu 65 (Mi = -5/2, -3/2, -1/2, 1/2, 3/2, 5/2) present. If the two copper 
isotopes coexist, they couple with one another giving rise to poor resolution in the form of M-shape and W- 
shapes as observed at the beginning and end of the g// lines. In this case, the isotope effect is taken into account 
in search for the best line of fit using computer simulation. The effective values of g, A and AH pp corresponding 

to the simulated EPR lines of Cu 2+ in CAS at different temperatures were recorded. However, the effect of 
isotope is neglected by setting a=l where there are no M-and W- shapes present in the EPR spectra. 



Results and discussion 
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Cu 2+ :CAS 



at different 



Fig. 1 Observed (solid lines) and simulated (dashed lines) spectra of 
temperatures 

The EPR spectra of Cu 2+ in CAS (Fig. 1), in an ac plane, show three groups of four resolved hyperfine 
lines in such a way that one set appeared at the low magnetic field and the remaining two at high magnetic field. 
The high resolution of the spectra in the high magnetic field region at 15 K results to reduction in broadening 
due to spin lattice relaxation time. The solid lines represent the observed spectra while the dashed lines are the 
simulated spectra. The complete anisotropic spectra of the hyperfine lines in the temperature range (15<T<180 
K) has been attributed to phonon-induced reorientation rate or quantum tunnelling of the Cu 2+ ions among the 
three JT 



Table 1 Observed and simulated g, A and AH PP values 



T(K) 


gz 


gx 


gy 


g<lll> 


A Z (G) 


A X (G) 


A V (G) 


A<111>(G) 


AH DD (G) 


15 


2.415 


2.144 


2.094 


2.222 


104 


20 


17 


62 


12.0 




2.416 


2.136 


2.093 


2.214 


102 


20 


16 


62 


12.0 


80 


2.401 


2.150 


2.101 


2.221 


96 


12 




56 


12.0 




2.397 


2.148 


2.101 


2.220 


95 


** 


** 


** 


12.0 


100 


2.393 


2.151 


2.103 


2.219 


94 




* 


* 


14.0 




2.393 


2.149 


2.105 


2.219 


93 


** 


** 


** 


14.2 


120 


2.384 








91 








16.0 




2.380 


** 


** 


** 


92 








16.3 


140 


2.370 








88 






* 


20.0 




2.367 


** 


** 


** 


88 








19.8 


180 


2.353 








80 








24.0 




2.348 


** 


** 


** 


80 








24.8 



where bold face values represent simulated values, * and ** represent missing observed and simulated values 
respectively. 
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(b) 

Fig. 2 Variation of observed and simulated (a) g-values and (b) A- values with temperature 
potential wells and Cu 2+ ions would then stabilize statistically among the three JT wells (Yerima and De, 2014). 
The observed and simulated values of the hyperfine coupling constant (A), g-factor (g) and peak-to-peak 
linewidth (AH PP ) at various temperatures are recorded in Table 1. The results show that A z and g z measured in 

the parallel direction of the magnetic field decrease with temperature while g±(g x ,g y ) and A ± (A x ,A y ) increase 
with temperature. 

The main difficulty in computer assisted structure elucidation consists in the extraction of the relevant 
spectroscopic parameters such as resonance frequency, coupling constants and linewidths from the spectra. In 
this study, this problem has not been solved in totality because of highly complex or insufficiently resolved 

splitting of the g± lines of Cu 2+ :CAS EPR spectra at high temperatures. However, the complete resolution of the 
hyperfine lines has provided reliable information that the number of equivalent nuclei in this system is three 
resulting to the four hyperfine lines of intensity ratios 1:3:3:1 in each of the three groups of four hyperfine lines. 



Free ion ground 
Term 



Crystal field Vibronic 
splitting coupling 



Strain's 
splitting 



Zeeman 



splittimg 




Fig. 3 Various splitting processes of 2 D ground state of cu 2+ 

Fig. 3 shows the representation of the effects of the cubic crystal field, vibronic coupling, strain 
splitting, and Zeeman interaction on a 2 D electronic state. Allowed AM S =±1, AM^O EPR transitions are 
indicated by straight arrows at the right hand side of the Figure. Curved arrows correspond to vibronic 
relaxation effects which can produce an isotropic spectrum by averaging a portion of the strained -broadened 
resonance pattern. A host of authors (Reynolds and Boatner, 1975; Boatner and Reynolds, 1977) have shown 
that for pure static JT systems like Cu 2+ :CAS, the ratio of the random strain splitting (tf) to the tunnel splitting 

(3r) is greater than 5 i.e. S > 5 X 31" . Thus, for the residual linewidth, AH PP =12 G and g=2.415 at 15 K, 
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5 > 5 x 3r = 5 x i^iif p p = 155 .88 G = 0.0176 cm" 1 and similarly at 180 K, AH PP =24 G and g=2.353, we 

have S = 0.0342 cm' 1 . Thus S increases with temperature and is far less than the average crystal splitting 
A=289.6 cm" 1 (Yerima and De, 2014) computed from intensity and linewidth measurements. 
In the case of Cu 2+ :CAS, the presence of M- and W-shapes at the beginning and end of the g// lines at low 
temperatures is attributed to the coupling of copper isotopes (copper 63(69%) and 65(31%)) (Yerima, 2007). 
The copper 63 ion has 3 equivalent protons while copper 65 ion has 4 which are expected to give 4 and 5 

respectively resolved g± or g// lines. At low temperatures, say 15 K, coupling of copper 65 ions becomes 
significant interacting with that of copper 63 ions giving rise to the M- and W-shapes. At high temperatures, 

copper 63 ions dominates giving rise to the four resolved g /f lines and the g ± lines are insufficiently resolved 
(Fig. 1). 

CONCLUSION 

In this paper, the theoretical expression for the line shape function of the orthorhombic JT spectra of 
Cu 2+ in CAS system have been reported. The spectra of Cu 2+ :CAS were simulated by the aid of computer using 
spectrum fitting technique at various temperatures. The results show that the simulated spectra match excellently 
well with the observed spectra especially at the low field region where the hyperfine lines are sufficiently 
resolved. The splitting processes in the 3d 9 orbital states of Cu 2+ have been discussed. The criterion of pure 

static JT systems demands that the condition — > 5 must be satisfied, where T = — - — and 8 = the 

random strain splitting was found to be about 0.0176 cm" 1 at 15 K and 0.0342 cm" 1 at 180 K. The insufficient 
resolution of the spectra especially at the high magnetic field side and temperatures is attributed to the rate of 
coupling of copper ions in this system. 
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